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EVALUATION 


The  purpose  of  this  effort  was  to  demonstrate  the  concept  of  Equivalent 
Aperture  Loading  (EAL)  and  interface  this  concept  with  the  "Method  of  Moments" 
(MOM)  analysis  called  "Bodies  of  Revolution"  (BOR).  The  reason  that  this 
combination  was  chosen  is  that  previously  MOM  and  BOR  were  applied  to  a 
conceptual  aperture  coupling  problem  on  board  a  missile.  The  results  were  found 
to  be  questionable  because  of  the  inability  to  accurately  model  the  complex  shaped 
aperture  involved.  This  effort  combined  the  EAL  technique  with  MOM  in  an 
attempt  to  correct  this  deficiency  and  at  the  same  time  demonstrate  that  such  a 
hybridization  of  techniques  is  cost  effective  in  terms  of  computer  resources. 

The  results  of  the  field  mapping  inside  the  missiie  are  contained  in  this  report 
and  will  be  compared  to  those  obtained  from  measured  data  associated  with 
another  RADC/RBCT  project  when  it  is  completed. 

The  advantage  of  this  technology  are  a  computer  cost  savings  by  hybridizing 
two  or  more  analysis  techniques  while  maintaining  a  reasonable  degree  of 
accuracy. 

DANIEL  E.  WARREN 
Project  Engineer 
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SECTION  1 


INTRODUCTION 

Bodies  of  revolution  (BORs)  are  rotationally  symmetric  structures  that  approxi¬ 
mate  the  enclosures  of  many  practical  systems  such  as  missiles  and  aircraft  fuselages. 
The  aperture  loading  technique  f ALT )  is  applied  here  to  predict  the  plane  wave  ex¬ 
cited  E~  and  H-fields  penetrating  a  thin,  conducting  walled  BOR.  The  BOR  has  an 
open  front  end  and  a  thin,  dielectric-filled,  sleeve-fit  seam  formed  by  overlapping 
walls.  The  present  version  of  the  computer  code  implementing  ALT  applies  to  arbi¬ 
trarily  shaped  BORs  with  any  number  of  simple  apertures  (circumferential  W3i1  sepa¬ 
rations)  and  one  thin  capacitive  seam.  The  code  can  apply,  with  minor  mod i f  teat  ion, 
to  any  number  of  odd-shaped  seams  and  openings  with  material  fill;  e.g.,  conductive 
putty,  fiber  glass,  or  composite. 

The  aperture  characteristics  are  dealt  with  naturally  in  a  localized  manner. 

In  the  case  of  composites,  this  permits  relative  ease  in  accounting  for  the.'  an  iso¬ 
trophy  of  the  medium.  The  aperture  characterise i"s  can  be  varied  for  reasonably 
small  apertures  with  only  minor  changes  in  computation.  This  is  particularly  use¬ 
ful  for  apertures  that  are  iii-dafined  where  only  an  average  respo.it-  or  other  sta¬ 
tistical  measure  is  appropriate.  An  ill-defined  aperture,  for  example,  is  a  g.Uo'  a- 
fit  seam  in  a  missile  shell;  the  spacing  between  seam  wails  may  vary  from  missile  to 
missile  since  the  fit  is  designed  to  be  as  "tight"  as  possible.  St.-ial  Samples 
such  apertures  are  shown  in  Figure  1. 

The  computer  code  extends  the  capability  of  an  earlier  BOR  apart ur-  o  ',  '.  ;- 

code,  B0R3  [9],  principally  by  employing  ALT.  Additional  extensions  and  i;  .ce¬ 
ments  include  near  H-field  computation,  P-  as  well  as  0-polarized  plane  wave  ex¬ 
citation,  computer  memory  requirement  reduction  by  approximately  60  percent,  u.d 
out-of-core  matrix  computation  and  inversion.  The  out-of-core  capability  ows 
large  problems  to  be  solved  if  long  run  times  can  be  tolerated.  The  present  version 
of  the  code  [4]  is  a  collection  of  several  codes  and  is  operational  on  a  PDP-11 ! 34  mini¬ 
computer.  The  solutions  presented  in  this  report  were  obtained  with  the  PDP-11 /54. 

As  with  B0R3,  the  ALT  BOR  code  is  based  on  a  Harr ington/Mautz  BOR  moment  method 
for  obtaining  plane  wave  excited  scattering  currents.  The  Harr ington/Mautz  code 
employed  in  developing  the  ALT  BOR  code  is  a  revision  of  the  original  Harrington/ 

Mautz  code  [12]  employed  in  B0R3.  The  revised  Harr ington/Mautz  code  [1]  sacrifices 
some  accuracy  in  matrix  element  computation  but  achieves  at  least  50  percent  faster 
run  time.  The  revision,  however,  is  applicable  only  to  closed  bodies;  therefore  an 
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SLEEVE  FIT  SEAM 


Figure  1.  Body  cf  Revolution  with  Several  Types  of  Apertures 


edge  element  expan function  Wda  incorporated  in  the  ALT  BOR  _ode  in  the  a  acne  wan¬ 
ner  as  in  the  original  "a  mngton /Mant  z  code. 

It  also  was  necessary  to  alter  the  near  E-field  computation  from  that  of  B0R3; 
however,  it  still  follows  closely  Bevensee’s  approach  [6]  as  employed  in  B0R3.  A 
near  H-field  computation  also  was  included.  The  theory  and  coding  were  adapted  from 
a  Harrington/Mauti'  H-field  formulation  BOR  moment  method  fll. 

A  third  Harrington/Muutz  E-field  formulation  moment  method  BOR  scattering  code 
has  recencly  become  available  [13].  This  code  is  superior  to  the  previous  Harrington/ 
Mautz  codes,  especially  in  the  manner  in  which  it  treats  edge  effects.  Following 
the  work  of  Wilton  and  Clission  [10],  they  employed  a  "shifted"  pulse  type  expansion 
function  for  the  i-directed  surface  current  component  with  excellent  success  in 
removing  the  computationally  induced  surface  current  oscillations  that  occur  with 
some  geometries.  An  ALT  BOR  code  based  on  this  improved  BOR  code  presently  is  under 
development . 

The  ALT  BOR  code  also  is  based  in  part  on  the  equivalent  aperture  excitation 
technique  (HAET).  The  EAET  is  characterized  by  greater  sensitivity  over  the  conven¬ 
tions  ;  -ear-field  scattering  method  [A]  in  that  small  field  strengths  are  determines 
more  accurately.  In  the  conventional  method,  the  surface  current  induced  on  the  con¬ 
ducting  portion  of  the  body  and  its  corresponding  scattered  field  are  determined. 

The  internal  field  then  is  found  by  adding  the  impressed  (exciting  field  in  the  ab¬ 
sence  of  the  body)  and  scattered  fields.  For  small  field  levels,  such  as  those  as¬ 
sociated  with  tiny  apertures,  this  computation  involves  a  subtraction  of  almost 
equal  numbers  with  a  corresponding  loss  in  accuracy.  This  difficulty  is  avoiccd 
with  the  EAET  because  the  internal  field  is  computed  as  a  radiation  problem  •■>nce  an 
equivalent  aperture  excitation  is  determined.  The  improvement  in  sensitivit-  of  ib" 
EAET  over  the  conventional  method  has  been  demonstrated  [9]. 

The  B0R3  code  was  the  first  to  employ  the  EAET.  The  EAET  subsequent Iv  lias  been 
employed  in  a  moment  method,  body  of  translation,  aperture  coupling  code  [14].  The 
EAET  also  has  been  incorporated  in  a  moment  method  BOR  code  in  a  more  efficient  man¬ 
ner  than  in  BOR3  [10l» 

A  detailed  theoretical  treatment  of  the  ALT  BOR  code  from  basic  BOR  theory  to 
the  EAET  is  given  in  Appendix  A.  ALT  is  a  natural  extension  of  the  EAET  and  the 
principal  reason  for  developing  the  ALT  BOR  code.  The  theoretical  basis  of  ALT  is 
given  in  Section  2. 
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SECTION  2 


PLANE  WAVE  EXCITATION  OF  BODIES  OF  REVOLUTION  WITH  COMPLEX  APERTURES 

An  E-field  formulation,  moment  method  solution  of  body  of  revolution  (BOR)  ap¬ 
erture  coupling  problems  is  detailed  in  Appendix  A,  Consider  a  plane  containing  a 
straight  line  and  an  arbitrarily  curved  line.  The  surface  of  a  BOR  is  traced  by  the 
curved  line  when  the  plane  is  rotated  with  the  straight  line  as  its  axis.  The 
straight  line  is  called  the  axis  of  the  BOR  and  is  typically  the  z  axis  of  a  standard 
p$z  cylindrical  coordinate  system.  The  curved  line  is  called  the  generating  curve 
of  the  BOR.  These  and  other  characteristics  of  a  BOR  are  discussed  in  Appendix  A, 
which  also  contains  a  description  of  the  equivalent  aperture  excitation  technique 
(EAET) . 

The  B0R3  computer  code  [9]  implements  EAET,  but  is  hampered  by  an  unnecessarily 
large  computer  memory  requirement;  in  addition,  B0R3  applies  only  to  9-polarized 
incident  waves  and  does  not  determine  internal  H-fields.  The  ALT  BOR  code  based 
on  the  theory  in  this  report  extends  B0R3  to  include  <j>-polarized  incident  waves 
and  near  H-field  computation.  It  also  requires  about  60  percent  less  computer  mem¬ 
ory  than  does  B0R3.  In  addition,  the  ALT  BOR  code  employs  the  aperture  loading  tech¬ 
nique  (ALT)  suggested  in  [9(Section  5.2)},  but  not  previously  implemented.  ALT  ap¬ 
plies  to  apertures  considerably  more  complex  than  simple  circumferential  openings. 

The  extension  of  EAET  to  include  ALT  is  discussed  below. 

Following  the  BOR  theory  of  Appendix  A,  the  E-field  component  tangential  to  ap¬ 
erture  surfaces,  EA  ,  is  expanded  as 
tan’  r 


EA 

'  tan 


E(EAfc  t  +  EA**1  <£) 
np  np 


np 


where 


:ac  =  -VA1  <5  (t-t  )  — - 
np  np  p  2  it 


jn* 


jn<t> 


EAr  =  -VA*  6 (t-t  )  ~ 
np  np  p  2r 


(1) 

(2) 

(3) 


The  curvilinear  coordinate  t  in  (l)-(3)  (Figure  A-l)  traverses  the  BOR  generating 
curve  (latitudinal  variation),  the  4>  coordinate  (Figure  A-l)  is  the  azimuthal 
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(longitudinal)  coordinate,  6 ( t )  is  the  unit  impulse  function,  the  integer n(-°® < n <«) 
spans  the  azimuthal  BOR  (sinusoidal)  modes,  the  interger  p  (p>0)  spans  the  BOR  gener¬ 
ating  curve  triangle  functions  that  reside  in  an  aperture,  tp  is  the  t  coordinate  of 
the  peak  of  the  p^  aperture  triangle  function,  and  VA^j  >  are  coefficients 

in  units  of  volts. 

The  surface  current,  JA,  flowing  through  the  apertures  when  the  apertures  are 
excited  by  (1)  is  expanded  as 

JA  =  y'  (IA*1  ft(t)ejn<f>  t  +  IA4*  f ^(t)e^n<^  4>)  (4) 

L-i  nP  P  np  p 

np 


where  f  (t)  and  f^(t)  are,  respectively,  the  t-directed  and  <}>-directed  unit  triangle 

functions  divided  by  the  BOR  radial  coordinate,  p  (p  is  understood  to  be  a  function 

of  t),  and  IA11  ,  IA^  are  coefficients  in  units  of  amperes, 
np  np 

Let  the  pth  elements  of  the  matrix  column  vectors  IA*1 ,  1A^,  VA*" ,  and  VA^  be 
v  n’  n’  n’  n 

IA^,  IA^  ,  VAnp’  an<^  VAn  ’  resPectively-  As  discussed  in  the  appendix  (Section  A-6) , 

tli  tVi 

the  n  mode  of  aperture  surface  current  coefficients  is  related  to  the  n  mode  of 

aperture  E-field  excitation  coefficients  by  the  matrix  equation 


K') 


The  matrix  containing  the  [YA^C]  etc.  submatrices  is  called  the  n^  mode  "aperture 
admittance  matrix."  The  elements  of  this  matrix  are  obtained  from  the  BOR  moment 
method  theory  detailed  in  Appendix  A. 


The  equivalent  aperture  excitation  technique  (EAET)  employs  (5)  in  obtaining 
aperture  port  voltage  source  vectors,  VA^  and  VA^,  from  "equivalent  aperture  port 

current  excitation  vectors"  iAfc  and  IA^.  The  IAC  and  IA^  are  known  (Equation  (A-122)) 

n  n  n  n 

from  the  surface  current  induced  on  the  BOR  after  the  apertures  are  covered  with  per¬ 
fect  conductors.  This  is  discussed  in  Section  A-6.  It  also  is  shown  in  Section  A-6 
that  excitation  of  the  BOR  with  an  equivalent  aperture  current  excitation  in  terms 
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of  IA^  and  IA*  or,  more  conveniently,  the  aperture  tangential  E-field  in  terms  of 

VA^  and  VA^  (the  apertures  are  replaced  and  the  external  exciting  field  removed) ,  re¬ 
sults  in  the  internal  fields  of  the  original  problem. 

If,  however,  the  aperture  has  complex  edges  (e.g.,  a  sleeve-fit  seam)  or  is  ma¬ 
terial  filled  (e.g.,  with  a  dielectric)  or  both,  some  current  will  flow  through  this 
aperture  "load"  when  excited  by  the  equivalent  current  source.  This  is  accounted  for 
with  ALT  by  first  identifying  the  relationship  between  EAtan  and  the  surface  current, 
JL,  flowing  through  the  load.  Corresponding  to  (4) ,  the  expansion  for  JL  is 


JL  =  V*  (IL*  f*(t)  ejn*  t  +  IL*  f*(t)  ejn<()  $) 
/  i  np  p  np  p 


Since  the  walls  of  the  BOR  are  thin,  EA  fills  the  aperture  load  surface  and  is 

~  tan 

related  to  JL  through  local  phenomena  relations;  e.g.,  [11] 


JL  =  (o  +  jw(e-£  ))  EA 


0"  -tan 


if  the  aperture  is  filled  with  a  material  of  finite  conductivity  (o)  and  permit¬ 
tivity  (e). 

Equations  (1),  (2),  (3),  and  (6)  and  relationships  such  as  (7)  yield 


[YL^]  [0]  VA* 

n  n 


where  IL*1  and  IL^  are  matrix  column  vectors,  the  pC^  elements  of  which  are  IL*"  and 
n  n  np 

IL^  ,  respectively,  [0]  is  the  null  matrix,  and  [YLCt ]  and  [YL^]  are  diagonal  matrices, 
np  n  n 

The  matrix  containing  the  [YL^C]  etc.  submatrices  is  called  the  aperture  load  matrix. 
Continuity  of  current  constrains  the  total  aperture  current,  JA,  to  be 


JA  =  -JS  +  JL 
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where  -JS  is  the  equivalent  aperture  excitation  current  given  by  (A-120) .  Equations 
(4),  (A-120),  and  (6),  combined  via  (9),  yield 


r  &  i 

n 

->  f 

IS 

n 

+ 

’ll* ' 

n 

i 

-e-  C 

+< 

tm 

_ 1 

-e-  c 

tm 

_ 1 

->  a 

11/ 

n 

Equation  (10),  with  (5)  and  (8),  yields 

riYA«]  [YAC*] 

n  n 

[YA*C]  [YA**] 

L  n  n 

The  solution  of  (11)  for  VA^  and  VA^  provides  EA^ ^  via  (1)  through  (3).  As  shown  in 
Appendix  A,  Section  6,  EA^^  can  be  applied  as  an  excitation  resulting  in  the  inter¬ 
nal  fields  of  the  original  problem. 

If  the  qth  aperture  port  is  an  electrically  small,  dielectrically-filled  sleeve- 
fit  seam  (Figure  2),  the  corresponding  elements  of  (8)  lead  to 


[YLtt:] 
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[0] 


[03 


VA 


VA 


IS 


IS 


(11) 
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tt 

"nqq 


(12) 
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EXTERNAL  REGION 


(13) 


Figure  2.  Portion  of  BOR  Generating  Curve  Containing  a  Sleeve-Fit  Joint 
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Since  the  aperture  width,  h,  is  much  less  than  X,  VA^  may  be  assumed  zero.  The 
^-directed  "load"  current,  IL^,  can  easily  "skirt"  the  thin  aperture  as  if  the  ap¬ 
erture  were  absent  and  cancel  the  equivalent  port  current  excitation,  -IS^1  ;  thus, 
by  (10) ,  IA^  »  0  and  (11)  is  reduced  by  one  equation  and  one  unknown  by  removing  the 

qth  rows  of  [YA^],  [YA^]»  and  [YL^],  and  the  qC^  columns  of  [YA1"^],  [Y$^  ],  and 
n  n  n  n  n 


[YI^] ,  and  the  qC^  elements  of  and  fs^. 


The  aperture  load  admittance,  YL  ,  for  the  sleeve-fit  seam  is  approximated  as 
th  nqq 

follows.  The  j  triangle  function  is  located  with  peak  at  the  gap.  From  (2),  the 
voltage  across  the  gap  V  ^  is  approximately 


/  < 


t  dt  =  -  VAC 
nq  nq  2ir 


The  surface  current  flowing  through  the  gap,  containing  both  displacement  and  polar¬ 
ization  current,  is  given,  from  (6),  by 


I  -  II*  fC (t  )  e3”* 
gap  nq  q  q 


or,  since  fC(t  )  =  —  (note  discussion  following  (4)), 

q  q  P 


t  ejn<j> 

I  =  IL  - - 

gap  nq  p 


The  gap  susceptance  per  unit  length  is 


"  jwCg 

gap 


where  the  gap  capacitance  per  unit  length  C  is  approximately 

S 


r  -  SL 

g  h 
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SECTION  3 


APPLICATION  TO  A  MULTIPLE  APERTURE,  PARTITIONED  CAVITY  BOR 

The  aperture  loading  technique  (ALT)  was  applied  to  a  BOR  defined  by  the  gener¬ 
ating  curve  of  Figure  3.  The  excitation  was  the  nose-on  x  polarized  unit  plane  wave 

E6  =  -e"jkz  i  V/m  (21) 

defined  with  reference  to  the  coordinate  system  shown  in  Figures  A-l  and  A-2  and 
where  k  is  the  wave  number  corresponding  to  300MHz.  The  body  had  two  apertures: 
open  front  end  and  a  sleeve-fit,  dielectric-filled  seam.  A  conducting  partition 
was  located  at  z  =  0.255m. 

The  partition  was  assumed  to  meet  the  outer  wall  at  z  =  0.255m  since  h  was  elec¬ 
trically  very  small.  This  formed  a  three-surface  junction  which  could  not  be  modeled 
with  the  ALT  BOR  code  without  significant  code  modification.  Variations  to  the 
Figure  3  BOR  multiple  surface  junctions  are  shown  in  Figures  4  and  5.  An  accurate 
solution  to  the  Figure  3  aperture  coupling  problem  was  achieved  by  properly  defining 
the  Figures  4  and  5  problems  and  combining  their  solutions  within  the  equivalent 
aperture  excitation  technique  (EAET)  framework  described  below. 

The  EAET  naturally  divides  into  three  parts.  Part  1  uses  BOR  scattering  theory 
(Sections  A-l  through  A-3)  to  obtain  the  external  scattering  current  with  the  aper¬ 
tures  removed  (shorted).  Here  the  Figure  4  BOR  (solid  line)  applies  exactly.  Part 
2  uses  BOR  theory  (Section  A-6)  to  obtain  an  aperture  voltage  excitation  by  computing 
an  aperture  admittance  matrix  and  combining  it  with  an  equivalent  aperture  current 
excitation.  The  current  excitation  is  the  negative  of  the  external  scattering  cur¬ 
rent  induced  on  the  aperture  shorts  in  Part  1.  (The  ALT  improves  this  aperture  volt¬ 
age  excitation  computation  by  adding  an  aperture  load  admittance  to  the  aperture  ad¬ 
mittance  matrix  as  described  in  Section  2.)  The  aperture  voltage  excitation  is  ob¬ 
tained  by  exciting  both  external  and  internal  regions  simultaneously.  Neither  Fig¬ 
ures  4  nor  5  modex  the  Figure  3  body  exactly  for  this  computation.  Since  the  body  is 
electrically  thin  at  300  MHz,  the  internal  fields  decay  rapidly  away  from  the  aper¬ 
ture  and  the  Figure  4  model  is  expected  to  be  adequate  for  computing  the  aper¬ 
ture  voltage  excitation. 

Part  3  of  the  EAET  uses  BOR  theory  (Sections  A-4  and  A-5)  to  compute  the  inter¬ 
nal  fields  by  applying  tho  aperture  voltage  excitation  to  the  BOR.  This  computation 
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results  in  the  same  internal  fields  when  applied  to  the  Figure  5  body  as  for  the 
Figure  3  body. 

This  scheme  of  using  the  Figure  4  body  for  Parts  1  and  2  of  the  EAET  and  the 
Figure  5  body  for  Part  3  was  expected  to  result  in  a  good  approximation  to  the  inter¬ 
nal  field  of  the  Figure  3  body  when  excited  by  (21). 

The  numbers  of  triangle  functions  on  the  Figures  4  and  5  generating  curves  were 
chosen  to  provide  the  desired  resolution  in  near  field.  This  was  critical  only  for 
the  Figure  5  geometry  since,  as  discussed  above,  the  Figure  4  geometry  was  not  used 
for  near-field  computation.  The  Figure  4  expansion  functions  must  be  sufficient  to 
resolve  surface  current  and  aperture  voltage. 

Only  the  n  =  1  BOR  mode  fields  and  currents  are  excited  by  axially  incident 

0 

plane  waves.  The  incident  field  E  given  by  (21)  is  6  directed  in  the  4>  =  0  plane. 
From  (A-55) ,  therefore,  the  induced  surface  current  has  the  form 


surface  current  due 
to  Ee  of  (21) 

where  t  is  the  generating  curve  coordinate  (Figure  A-l) .  Surface  currents  computed 
by  the  ALT  BOR  code  are  prone  to  exhibit  severe  oscillation  particularly  if  the  BOR 
surface  has  edges  and  the  BOR  is  thin.  This  was  experienced  by  Davis  and  Mittra  [15]. 
The  problem  was  analyzed  by  Glisson  and  Wilton  [16],  and  they  developed  a  corrected 
BOR  moment  method  code.  Another  corrected  BOR  moment  method  code  recently  was  de¬ 
veloped  by  Mautz  and  Harrington  [13].  The  Mautz  and  Harrington  version,  more  closely 
resembling  the  BOR  code  upon  which  the  ALT  BOR  code  has  been  developed,  is  in  the  proc¬ 
ess  of  being  adapted  to  the  ALT  BOR  code. 

The  present  version  of  ALT  BOR  code,  meanwhile,  has  been  applied  to  the  Figure 
3  problem  via  the  Figures  4  and  5  geometries  as  described  above.  The  computed  scat¬ 
tering  current  on  the  Figure  4  geometry  exhibited  severe  oscillations  in  both  Jt  and 

A 

J  .  In  attempting  to  reduce  this  oscillation,  the  tail  section  of  the  Figure  4  geo¬ 
metry  was  rounded  as  indicated  by  the  dashed  line  in  Figure  4.  This  avoided  the  edge 
at  reference  point  7.  The  curved  path  was  chosen  such  that  the  path  lengths  6a  to  8a 
and  6a  to  7  to  8  were  equal.  The  generating  curve  then  was  segmented  (for  triangle 
expansion  functions)  in  a  manner  maintaining  less  than  2:1  variation  in  adjacent  seg¬ 
ment  sizes.  The  rounded  tail  section  and  small  variation  in  adjacent  segment  sizes 
resulted  in  the  J1"  and  J1**  curves  of  Figure  6.  The  generating  curve  was  divided  into 


J  cos  (}>t+Jysin(J><t> 


(22) 
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Continued 


86  segments  (42  triangle  functions).  Thc  J  current  then  was  reasonably  smooth.  The 
J '  current,  however,  still  oscillated  severely.  Nevertheless,  the  ALT/EAET  was  con¬ 
tinued  in  the  belief  that  the  internal  fields  were  not  strongly  affected  by  the  oscil 

f. 

lation  in  J'.  In  Figute  5,  the  currents  are  plotted  vs.  the  z  coordinate  of  the  gen¬ 
erating  curve.  Each  plot  point  corresponds  to  the  peak  of  a  triangle  function. 

Two  sets  of  five  aperture  voltages  sources  (four  at  the  front  end  and  one  at 

the  sleeve  fit  seam)  were  determined.  One  set  corresponded  to  the  air  filled  seam 

and  one  to  the  dielectric  filled  seam  (dielectric  constant  c  =  5.5).  These  sources 

r 

were  applied  to  the  Figure  5  generating  curve  as  indicated  in  Figure  5,  and  the  inter 
nal  E-  and  H-fields  computed.  The  Figure  5  generating  curve  was  divided  into  82  seg¬ 
ments  (40  triangle  functions). 

The  internal  fields  have  the  form  (note  (A-78)  and  (A-118)) 

E-field  due  to  , 

e  ,  ....  ?  =  E  cos  0  p  +  E  cos  0  z  +  E,  sin  f  $  (23) 

E  of  (21)  >  p  z  <f> 

H- fie  Id  due  to  . 

_e  r  ....  ^  =  H  sin  <|>  p  +  H  sin  ;  z  +  H  cos  ■+•  $  (24) 

E  of  (21)  f  p  z  } 

Contour  maps  of  the  magnitude  of  E  ,  E  ,  E_,  r.H  ,  -H  ,  and  ■!  in  dB  are  shown  in 

P  Z  fi  i  z 

Figures  7  and  8  where  r,  is  the  free  space  wave  impedance.  The  resolution  of  computed 
internal  field  data  to  achieve  these  plots  was  0.005  m  within  0.02  m  of  an  aperture 
and  0.01  m  elsewhere.  The  field  contours  were  smoothed  wherever  oscillations  were 
likely  caused  by  the  granularity  of  the  surface  current  representation.  Field  oscil¬ 
lations  were  evident  primarily  within  a  centimeter  or  so  of  a  conducting  surface. 

The  fields  (Figure  8)  in  the  dielectric  filled  seam  case  were  approximate¬ 
ly  25  dB  below  those  (Figure  9)  for  the  air  filled  seam  case  near  the  seam. 

This  could  have  been  predicted  by  comparing  the  n  =  1  mode  radiation  admittance 
(YA^*",.  j.  of  (11))*  at  the  seam  with  the  aperture  load  admittances  (YL^t^  of  (20)) 
for  the  air  filled  seam  and  the  dielectric  filled  seam  as  follows.  The  aperture  vol¬ 
tage  sources  were  determined  from  the  equivalent  aperture  current  sources  (the  latter 
are  independent  of  aperture  loading)  by  adding  together  the  radiation  and  load  admit¬ 
tances  and  inverting  the  sum  as  suggested  by  (11).  In  the  problem  under  consideration 

*  The  5,5  matrix  element  corresponds  to  the  seam  aperture  for  this  problem. 
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YA“  ,  -  0.04  -  jl.34  mS 

and  for  the  air  filled  seam  and  dielectric  filled  seam  respectively 

<5,5  =  j1-78*5 

r  =  1.0 

=  j  9 . 82  mS 

r  =  5.5 

Thus,  the  total  admittance  in  each  case  was 


air 


YA 


tt 

1,5,5 


•k,.{ 


=  0.04  +  j  0.44  mS 


1.0 


diel 


=  YA 


tt 

1,5,5 


+  <5,5[ 


5.5 


0.04  +  j  8.48  mS 


The  smaller  Y^  resulted  in  a  correspondingly  larger  aperture  voltage  excitation. 
The  air  filled  seam  nearly  "resonated"  the  aperture  by  almost  cancelling  the  radia¬ 
tion  inductive  reactance.  The  dielectric  filled  seam,  on  the  other  hand,  introduced 
a  shorting  effect  by  adding  a  large  aperture  capacitance  to  the  relatively  small 
aperture  radiation  reactance.  Since 


20  log 


Y  /  Y  ' 
1  air'  diel 1 


-25.6  dB 


the  '25  dB  smaller  fields  for  the  dielectric  filled  seam  than  for  the  air  filled 
seam  was  to  be  expected. 

The  air  filled  seam  problem  has  been  analyzed  previously  by  a  time  stepping  solu¬ 
tion  to  the  time  dependent  differential  Maxwell  equations  [17].  Maxwell’s  equations 
were  solved  at  each  of  many  time  steps  until  a  steady  state  solution  resulted.  Their 
results  differ  significantly  from  ours.  Some  of  this  difference,  perhaps,  can  be 
attributed  to  the  unwarrented  oscillations  in  surface  current  as  computed  with  the 
present  ALT  BOR  code.  The  ALT  BOR  code  is  presently  undergoing  significant  improve¬ 
ment  in  this  area  as  mentioned  above. 
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The  ALT  BOR  code  was  implemented  on  a  DEC  PDP-11/34  minicomputer.  The  above 
problem  required  hours  of  central  processing  time  with  this  machine  due  to  the 
enormous  amount  of  input/output  time  necessary  to  accomodate  large  complex  matrices 
(orders  approximating  80)  with  128K  bytes  of  main  memory. 

The  equivalent  approximate  run  time  assuming  the  entire  problem  was  solved  with 
main  memory  (no  peripheral  storage)  is  given  by 

333  222  333 

Time  =  n(A(m^  +  m^  +  a  )  +  B(m^  +  nt,  +  a  )  +  Ci(m^  +  m2  +  a  )  +  Dim^  +  Epn^) 

(25) 


where 


m^  =  number  of  triangle  functions  on  the  exterior  problem  generating  curve 
(e.g.,  that  of  Figure  4) 

m2  =  number  of  triangle  functions  on  the  interior  problem  generating  curve 
‘  (e.g.,  that  of  Figure  5) 

a  =  number  of  triangle  functions  within  the  apertures 

p  =  number  of  near  field  points 

i  =  number  of  different  incident  plane  wave  directions 
n  =  number  of  BOR  modes 

A,  B,  C,  D,  and  E  are  system  dependent  constants. 

A  second  computation  with  a  change  in  aperture  loading  (e.g.,  adding  dielectric  ma¬ 
terial  to  the  sleeve-fit  seam  of  Figure  3)  requires  only  the  additional  time  corres- 

3  2  2 

ponding  to  doubling  Aa,  Ci(i»2  +  a  ),  and  Epi^ . 

The  coefficients  in  (25)  for  a  Honeywell  6000  series  system  are 
A  =  3.2xl0~3 
B  =  1 . 3xl0~2 
C  «  5.0xl0_5 
D  =  3.2xl0"2 
E  =  0.87 

where  the  near  field  computation  time  coefficient  E  is  for  total  E-field  (scattered 
plus  impressed) . 

2  2 

The  principal  computer  main  memory  constraint  is  m^  +  m2  complex  words,  or  if 

only  the  "exterior  or  "interior"  problem  matrix  is  in  main  memory  at  one  time,  the 
2  2 

larger  of  m^  and  n^. 
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APPENDIX  A 


BODY  OF  REVOLUTION  SCATTERING  AND  APERTURE  COUPLING 


This  appendix  contains  the  description  of  an  extension  to  the  Mautz/ 

Harrington  [1]  moment  method  solution  of  an  E-field  formulation  for  scattering  from 
closed  conducting  bodies  of  revolution  (BORs).  The  following  conditions  are 
applicable: 

1.  The  BOR  may  have  edges  (disks,  tubular  cylinders,  etc.), 

2.  The  BOR  may  have  apertures  that  are  rotationally  symmetric  about  the 
BOR  axis , 

3.  An  electrically  thin  aperture  may  be  of  complex  shape  (e.g.,  "sleeve"  fit) 
or  material  filled, 

4.  The  scattered  near  electric  fields,  internal  and  external  to  the  BOR,  can 
be  computed,  and 

5.  The  scattered  near  magnetic  fields,  internal  and  external  to  the  BOR,  can 
be  computed. 

The  development  follows  closely  that  of  Mautz  and  Harrington. 

It  is  desired  to  determine  the  electric  surface  current  and  the  near  scattered 
field  of  a  perfectly  conducting  BOR  (Figure  A-l)  excited  by  an  incident  plane  wave. 

In  Figure  A-l,  p,  <J>,  and  z  are  cylindrical  coordinates  and  t  and  <)>  form  an  orthogonal 
curvilinear  coordinate  system  on  the  surface,  S,  of  the  BOR.  The  t  coordinate  tran- 
verses  the  "generating  curve"  of  the  BOR.  t  and  $  are  orthogonal  unit  vectors  in 
the  t  and  <j>  directions,  respectively.  The  coordinate  origin  is  on  the  axis  of  the 
BOR  but  not  necessarily  at  the  lower  pole  as  in  Figure  A-l.  Figure  A-2  defines  the 
propagation  vector,  k^,  of  the  incident  plane  wave  and  the  transmitter  coordinate  0g. 
Note  that  the  transmitter  <}>  coordinate  is  zero  so  that  kg  lies  parallel  to  the  xz 
plane.  In  Figure  2,  §  and  <}>  are  unit  vectors  in  the  0  and  y  directions,  respec- 
tively . 

Consider  separately  a  9  polarized  incident  plane  wave  defined  by 


-jk  '  T 

e  -  e 

E  =  6  e 
e 


and  a  polarized  incident  plane  wave  defined  by 


-jk  -r 

p  J  -  p  - 

E  =  e 
e 


(A-l) 


(A-2) 


A-l 


where  E6  denotes  the  incident  electric  field  and  r  is  the  radius  vector  from  the  ori¬ 
gin.  Either  plane  wave  gives  rise  to  t  and  <p  directed  electric  surface  currents 
on  S. 

1.  E-FIELD  SOLUTION 

The  E-field  solution  is  obtained  by  applying  the  method  of  moments  to  the  E-field 
integral  equation.  The  E-field  integral  equation  is  derived  by  setting  the  component 
tangential  to  S  of  the  total  electric  field  equal  to  zero  on  S. 

The  boundary  condition  that  the  total  tangential  electric  field  is  zero  on  S  is 
written  as 


-ES  =  Ee  on  S  (A-3) 

'tan  'tan 

S  6 

where  E  is  the  electric  field  due  to  the  electric  surface  current  on  S  and  E  is 
the  incident  electric  field  given  by  either  (A-l)  or  (A-2).  The  subscript  tan  de¬ 
notes  tangential  components  on  S. 

The  field  ES  can  be  expressed  in  terms  of  a  vector  potential,  A,  and  a  scalar 
potential,  $,  as 

ES  =  -  juiA(J)  -  V4>(J)  (A-4) 


where 

A(J) 


II  J(r')  -  ds' 

•'/  4  it  |r-r'  I 


®(J) 


i  ft.  *-iklr ti  ds. 

ZJJ  4-rr  |r-r’  I 


(A-5) 


(A-6) 


Here,  r  and  r’  are  vectors  to  the  field  and  source  points,  respectively,  J(r’)  is 
the  electric  surface  current  on  S,  k  is  the  propagation  constant,  u  is  the  permea¬ 
bility,  e  is  the  permittivity,  and  a  is  the  surface  charge  given  by 


a 


~  lim 
AS-H3 


^  J(r) -ndc  j 


~r~  V  -J(r) 
ju  'S 


(A- 7) 


A-3 


where  n  is  the  unit  vector  tangential  to  S  and  normal  to  the  curve,  C,  which  bounds 
the  small  portion,  AS,  of  S.  n  points  away  from  AS.  The  operator  V  •  is  the  sur¬ 
face  divergence  on  S. 

Following  the  method  of  moments,  J  is  expanded  as 


J  =  3^  (lt,  JC.  +  I*.  J*.  ) 
„  ,  \  nj  ~nj  nj  ~nj  / 
n  »  J 


(A-8) 


t  d>  t  <t> 

where  I  ,  and  I  .  are  coefficients  to  be  determined  and  J  .  and  J  .  are  expansion 
nj  nj  ~nj  ~nj 

functions  defined  by 


ii 


tfj(t)  eJn« 


(A-9) 


-nj 


$f*(t)  ejn<l 


The  choice  of  functions  f.(t)  and  f;(t)  is  discussed  below.  The  integral  over  S  of 

3  J  £ 

the  dot  product  of  (A-3)  with  each  one  of  a  collection  of  testing  functions,  W  , 

^mi’  by 


ii  -  £i5«>  ''lm* 


«.i  ■  e'Jm* 


(A-10) 


is  taken  to  obtain  the  matrix  equation 

?t  +  M  *♦)  =  vc 

a— 4  \  [_  mnj  n  |_  mnj  n  /  m 


(A-ll) 


E([c]  *:♦  Or]  r)  ■*, 


m 


th 


where  the  [ Z ] s  are  square  matrices  whose  ij  elements  are  defined  by 


ZP<1 

mnij 


f/< 1  -  (  +  ?»<£,>)  "S 


(A-12 ) 


A- A 


where  p  may  be  either  t  or  4>  and  q  may  be  either  t  or  if  .  Also  ,  A  and  A  are 

th  mm 

column  vectors  whose  i  elements  are  given  by 


-ffi 


,P  •  Eeds 
mi 


(A-13) 


where  p  may  be  either  t  or  <f>.  Lastly,  I  and  1^  are  column  vectors  of  the  coeffi- 
t  4  n  n 

cients  I  and  A  appearing  in  (A-8). 

The  following  manipulations  serve  to  transfer  the  differential  operator  on  4> 

in  (A-12)  to  WP . .  If  S  is  closed, 

-mi 


/A  ■ 


(<f>  W)  ds  =  0 


(A-14) 


where  W  denotes  W  If  S  has  an  edge,  the  surface  divergence  theorem  yields 


jj  •  ($  w)  ds  =  j  £• 


4>  W  dc 


where  C  is  the  path  of  the  edge.  Hence,  with  W  chosen  such  that  fW  •  0  along  the 
edge,  (A-14)  holds  whether  S  is  open  or  closed. 

The  representation 


v=  •  W  =  7TT  (p  W-  t)  +  ~  —  (W  •  ♦) 
-S  p  at  P  a  <p 


(A-15) 


of  the  surface  divergence  and  the  definition 


7  4>  =  t  +  4>  -~ 

'S  3t  p  3<(i 


(A-16) 


of  the  surface  gradient  imply  that 


V  •(<t’W)  =  $V-W  +  W-  V$ 

-g  ~  ~  ~  S 


(A-17) 


The  surface  gradient  of  in  (A-17)  can  be  replaced  by  the  ordinary  three-dimen¬ 
sional  gradient  of  <t>  because  (A-16)  is  the  component  of  the  three-dimensional  gradi¬ 
ent  tangential  to  S  and  W  is  tangential  to  S.  Substitution  of  (A-17)  into  (A-14) 
and  then  (A-14)  into  (A-12)  yields 


*pq-< 

mnij 


i  ■  -  °h  ♦<<&)) 


(A-18) 


where 


aP  =  -  V  •  WP 
mi  ju  -s  -mi 


(A-19) 


Since,  as  shown  in  Appendix  B  of  Reference  1,  (A-18)  is  zero  for  m  ^  n,  (A-ll) 
reduces  to 


w 

>  1 
n 

If 

>  ' 
n 

KJ 

'hi 

D  -©- 

1 _ 

i 

-e-  c 
f> 

_ i 

,  n=0 ,  ±1,  ±2, 


(A-20) 


where  ra  is  Jz^^j of  (A-18) .  It  also  is  shown  in  Appendix  B  of  Reference  1  that  the 
elements  of  Jzp<^j  are  given  by 


Z ^ . 

mi 


?+*■ 

"nij 


7^ 

nij 


"nij 


=  j  r\  j~  dt  J "  dt1  (  !<2pf*T(t)  p  ’f^(t')  (G^  sin  v  sin  v'  +  cos  v  cos  v') 

■  IE  (PfJ(t))  p’fj(f))  Ga)  (A-21) 

=  -n  f dt  pf^(t)  J dt'  (  k2p'fj (t')  G&  sin  v'  +  j  — —  {p'fj(t') )  G4  ) 


(A-22) 


=  n  y"dt  /dt'  P ' f ^ (t * )  (  k2pf^(t)  G&  sin  v  +  (  pf^(t) )  G^  )  (A-23) 


2 

jn  /  dt  pfj(t)  j  dt'p'fj(t')  (k2G5  —  ~~r  G^  ) 


(A-24) 


A-6 


I 


where  n  is  the  wave  impedance,  /p/c  ,  v  is  the  angle  between  the  tangent  to  the  gen¬ 
erating  curve  and  the  z  axis, 


■/ 


-j  kR 

d<f>'  ~kR~  cos  (n<f>') 


(A-25) 


S  = 


■/ 


it  -jkR 

d<J> 1  — cos  <t>'  cos  (n4> '  ) 


(A-26) 


f 


e~jkR 

d(j>'  — sin  *  sin  (n<$>') 


(A-27) 


and 


=  ^  (p-p')2  +  (z-z')2  +  4pp'  sin2(-^-  ) 


Here,  p,  z,  and  v  depend  on  t  while  p' ,  z',  and  v'  depend  on  t'. 

To  evaluate  (A-21)  through  (A-24),  pfP(t)  is  defined  as  a  four-impulse  approxi¬ 
mation  to  a  triangle  function  in  the  following  manner.  Letting  t  »  (p,z),  denote 
that  p  and  z  are  cylindrical  coordinates  of  the  point  t,  an  odd  number  greater  than 
or  equal  to  5  of  consecutive  points  t  =  (p,z)  *  t^  ■  (p^,  z^) ,  i  *  1,  2 ,  ...  P  on  the 

generating  curve  of  the  body  of  revolution  such  that  (p^,  z^)  and  (p~,  z~)  are  the 
poles.  If  the  body  of  revolution  has  no  poles  because  the  generating  curve  closes 
upon  itself,  as  with  a  torus,  then  two  points  must  be  overlapped  such  that 


(Pp_2 »  zp_  ♦  z2^  «ind  (pp>  zp)  (p2»  ^2^ 


The  generating  curve  is  approximated  by  drawing  straight  lines  between  the  points 
1,  2,  ...  P  and  defining  points 


( P 2 »  Z2>  >  i 


t  =  ^  =  (p^  z±)  = 


pi  +  pi+l  Zi  +  Zi+1 


(A-28) 


on  this  approximate  generating  curve.  The  length,  d^,  of  the  interval  centered 
about  t^  is  given  by 


A- 7 


t 


t2i-l  t2i-l 


fc21  t2i  t2i+l  ^2i+l  t2i+2t2i+2  t2i+3 


Figure  A-3.  Triangle  Function  (Solid)  and  Four  Impulse 
Approximation  (Arrows) 


(A-29) 


(pi+l  -  PI)2  +  (zI+ 1  '  z?2 


In  terms  of  coefficients  Tp+4i  4  defined  by 


kd 


2i-l 


4i-3  2<d2i.1+d21) 


41-2 


k(d21-l  +  2  d2l’d21 
d2i-l  +  d2i 


4i-l 


k(d2i+2  +  2  d2i+l)d2i+l 
d2i+l  +  d2i+2 


kd2  i+2 


4l  2(d2i+l  +  d2 i+2  ^ 


(A-30) 


one  constructs 


•'!<«  ■  i  S  1 

p=l 


p+4i-4  6(t_tp+2i-2) 


(A-31) 


where  6(t)  is  the  unit  impulse  function.  The  right-hand  side  of  (A-31)  is  the  de¬ 
sired  four- impulse  approximation  to  a  triangle  function  (Figure  A-3) . 

The  pff(t)  are  chosen  equal  to  pf^(t)  except  at  an  edge  where  the  4>  component 
of  surface  current  approaches  infinity;  therefore,  the  pf^(t)  functions  nearest 
edges  are  chosen  as  ramp  functions  (Figure  A-h).  A  four- impulse  approximation  re¬ 
sults  in 


PfJ(t) 


l  V 

k  h 


TRp+4i-4  ‘5(t-tp+2i-2) 


(A- 32) 


where  TRp+41_4  =  Tp+4i_4  except  for  TR 

of  the  generating  curve  or  for  T^ne_^ 
generating  curve  where  ne  =  number  of 
ginning  edge, 


^  and  TR2  if  an  edge  exists  at  the  beginning 

and  TR.  if  an  edge  exists  at  the  end  of  the 
4ne 

pfp  traversing  the  generating  curve.  For  a  be 


A-9 


0 


a.  Leading  Edge 


b.  Terminating  Edge 


Figure  A-4.  Triangle  Function  (Solid)  and  Four  Impulse 

Approximation  (Arrows)  for  4>  Directed  Current 


A-10 


(A- 33) 


Tf^  =  2kd1  -  T1 
TR2  =  2kd2  -  T2 
F  r  an  end  edge. 


TR,  =  2kd  .  -  T.  , 
4ne-l  ng-1  4ne-l 


2kd  -  T, 
ng  4ne 


(A- 34) 


where  ng  =  number  of  segments  traversing  the  generating  curve. 

For  ■—  (pf^(t)),  the  four  impulse  approximation 
dt  l 


£  (pfi(t)) 


4 

^  Tp+4i-4  6 (t-tp+2i-2) 


(A-35) 


to  the  derivative  of  the  triangle  function  of  Figure  A-3  is  chosen.  Figure  A  5  iT 
lustrates  A-35. 


1 


2i+1  +  d2i+2  | 

Figure  A-5.  Derivative  of  Triangle  Function  (Solid)  and  Four  Impulse 
Approximation  (Arrows) 


A- 11 


The  coefficients  T'p+4i_^  appearing  in  (A-35)  are  given  by 

T’  *  - ?.  -~J - 

41-3  d2i-l  +  d2i 


T!  .  =  — - - - 

41-2  d2i-l  +  d2i 


(A-36) 


41-1  d2i+l  +  d2i+2 


T,  -  2i+2 

41  ~  d2i+l  +  d2i+2 


where  d^  is  defined  by  (A-29) . 

Substituting  (A-31) ,  (A-32),  and  (A-35)  into  (A-21)  through  (A-24)  yields 
4  4 

Znij  =  Tn»  (Gs  sin  vi '  sin  V4*  +  ga  cos  v-t*  cos  v,i) 


-*4f(V  V  (G5  Sin  V  S±n  vj*  +  cos  r  cos  v ^ , ) 


p=l  q 


-  ^q’  G4> 


(A-37) 


4  4 

ZnSi  =  ^  X)  £  <TV  T  ,  (G,  sin  v . ,  +  r£-  TR.  T'  G.) 


p=l  q: 


_ j  \  »  *  I  v  .  ■  i  .  |  x  | 

pj  P  q  6  j'  kp^,  p'  q' 


(A-38) 


4  4 


TR  ,  G.  sin  v . ,  +  T\  TR  ,  G, ) 

q  6  i  kpj,  p'  q’  4 


(Ge  ~ 


k  o  .  ,o  . , 
i  y 


(A-39) 


(A-40) 


where 

p'  =  p  +  4i  -  4 
q’  -  q  +  4j  -  4 


(A-41) 


i'  =  p  +  2i  -  2 

j’  -  q  +  2j  -  2 


(A-41) 


The  subscript  i*  denotes  evaluation  at  t^,.  The  subscript  j’  denotes  evaluation  at 
t.,.  G.  ,  Gc,  and  G,  are  given  by  (A-25),  (A-26),  and  (A-27)  in  which  R  is  evaluated 

J  4  j  b 

at  t  =  t^,  ,  t'  =  tj ,  which,  in  terms  of  cylindrical  coordinates,  is  at  p,  z,  p' ,  z' *= 

P.,,  z.(,  p.,,  z.,.  If  i'  =  j 1 ,  R  is  replaced  by  an  equivalent  distance,  R  , 
i  l  J  3  ® 

given  by 


R  = 
e 


/ 


Cd.,/4)' 


2 

V 


(<j)’  /2) 


(A-42) 


(A-42)  may  be  obtained  by  displacing  the  field  point  a  distance  d^,/4  perpendicular 

to  the  plane  of  the  source  loop.  Now,  d.,/4  is  the  equivalent  radius  [2]  of  a  flat 

1  2 

strip  whose  width  is  d^, .  (A-42)  also  can  be  obtained  by  averaging  R  for  field 

points  displaced  a  distance  d^/4  in  either  direction  along  the  approximate  genera¬ 
ting  curve  from  the  source  loop. 


An  N,  point  Gaussian  quadrature  formula  is  used  to  calculate  the  integrals 
9 

G^,  G^ ,  and  G^  defined  by  (A-25),  (A-26),  and  (A-27).  According  to  this  quadrature 


formula, 

.IT 


L 


f(V)  d<ji' 


(A-43) 


where  f  ( <J> * )  is  the  function  being  integrated  and  x^  and  A^  are  constants  tabulated 
by  Krylov  [3].  In  (A-43),  the  multiplier,  j  ,  and  argument,  j  (xfc  +  1),  instead 


of  just  x^  are  due  to  the  transformation  of  Krylov's  interval  from  -1  to  1  into  the 
interval  from  0  to  tt. 


Since  replacement  of  i,  j,  p,  and  q  in  (A-41)  with  j,  i,  q,  and  p,  respectively, 
implies  replacement  of  i' ,  j',  p' ,  and  q'  in  (A-37)  through  (A-40)  with  j',  i',  and 
q' ,  and  p'^ respectively,  and  since  G^ ,  G5>  and  Gfe  are  symmetric  in  i’  and  j',  it  is 
evident  that 


(A-44) 


A-13 


(A-44) 


An  efficient  method  of  computing  (A- 34)  through  (A-40)  that  takes  advantage  of  (A-44) 
is  employed  in  the  computer  code  described  in  [4]. 

2.  PLANE  WAVE  EXCITATION 


th 

The  n  mode  generalized  voltage  vectors,  V  and  V  ,  are  matrix  column  vectors 

n  n 

with  elements  defined  by  (A-13)  (with  m  replaced  by  n) .  For  plane  wave  excitation, 
(A-13)  becomes,  after  substitution  from  (A-l) ,  (A-2),  and  (A-10) , 


Vnl  '/ dt  / 


_  „  -j  (k  -r+n<f>) 

dtp  (t-0  )  e 
e 


Vni  "  f dt  ofi^)  f  d*  (^§e) 


-j  (ke-r+n<(>) 


(A-45) 


.ydt  pfj(t>  f  d*  (t^e) 


-j  (k^r+nifO 


Z.  II 

V^i  =/dt  Dfi(t>  f  d<*>  (*4e) 


-j  (ke-r+nij>) 


where  the  second  superscript  on  V  is  0  for  the  incident  field  Ee  given  by  (A-l)  and 
4)  for  Ee  given  by  (A-2). 

With  a  view  toward  evaluation  of  (A-45),  it  can  be  seen  from  Figures  (A-2),  (A-6) , 
and  (A-7)  that 


t*0  =  -sin  0  cos  v  +  cos  9  sin  v  cos 

e  e  e 

<p’d  =  -cos  0  sin  <4 
e  e 

t*  9  =  sin  v  sin  <f> 

e 


A-46) 


A- 14 


(A-46) 


■  <p  =  cos 

e 


k  *r  ■  -  kz  cos  0  -  kp  sin  8  cos 

~e  '  e  e 


The  substitution  (A-46),  (A-31),  and  (A-32)  into  (A-45)  and  the  integral  formula 


,-n  r2v  j  (kp  sin  0  cos  <p  +  n<j>) 


J  (kp  sin  0  )  = 
n  e  z 


rf 


(A-47) 


deduced  from  (9.1.21)  of  [5]  for  Bessel  functions,  yields 


t0  irjn  \  *  T  ...  .  (-2J  sin  0  cos  v  +  j  (J  ,,  -  J  ,)cos  6  sin  v)e 
•  .  p+4i-4'  n  e  J  n+1  n-1  e 


jkz  cos  9 

J  c 


4^ 

V*®  =  T  .  .  .  (J  .  +  J  ^cos  0  e 

ni  k  *— i  p+4i-4  n+1  n-1  e 

P*1 


jkz  cos  9 


(A-48) 


t  <t>  irjn  \  ' 

V  ■  -  — /  .  T  .  /  •  /  (J  . ,  +  J  ,)sin  v  e 
ni  k  p+4i-4  n+1  n-1 

■  p=l 


jkz  cos  0 


rT  (j  -  j  )  e 

ni  k  p+4i-4^  n+1  n-1' 

p=l  v 


jkz  cos  B 


where 


J  =  J  (kp  sin  6  ) 
n  n  e 


(A-49) 


In(A-48),  p,  z,  and  v  are  to  be  evaluated  at  t  =  t  +2^_2‘ 


A-16 


3.  SCATTERING  CURRENT 


The  solution  to  (A-20)  yields  the  electric  surface  current  coefficients  IC  ,  I^1 . 

nj  nj 

for  each  BOR  mode  n.  The  expression  (A-8)  for  the  corresponding  surface  current  can 
be  simplified,  for  plane  wave  excitation,  by  combining  the  n  and  -n  terms.  The  sub¬ 
stitution  (A-9)  into  (A-8)  yields 


OQ 

jq  .  y  {(jt  >tq)£  +  (f4>  pq  ~} 

n  n 

n=-“> 


(A-50) 


where  fC  is  the  transpose  of  the  column  vector  of  the  f^(t),  f^  is  the  transpose 

of  the  column  vector  of  the  ff(t),  and  ItC^  and  1*^  are  column  vectors  of  the  co- 

1  n  n 

efficients  1^  and  I4^  respectively.  The  additional  superscript,  q,  is  either  8  or 
4>,  depending  on  whether  the  incident  electric  field  is  0  polarized  as  in  (A-l)  or  <)> 
polarized  as  in  (A-2).  The  column  vectors,  1^  and  1^ ,  are  obtained  by  solving  the 
matrix  equation  (A-20)  with  the  additional  superscript  q  on  the  column  vectors  there¬ 
in  to  denote  the  polarization  of  the  incident  electric  field. 

Inspection  of  (A-21)  through  (A-27)  reveals  that 
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It  is  apparent  from  (A-48)  that 
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Since  the  property  (A-51)  survives  matrix  inversion 
(A-52)  that 


n  =  0,  1,2,...  (A-51) 


n  =  0,  1,  2,  ...  (A-52) 


,  it  is  evident  from  (A-20)  and 


A-17 


n  *  0,  1 ,  2 ,  . . . 


(A-53) 


In  view  of  (A-53) ,  (A-50)  becomes 


(A-54) 


Je  *  (fc  I*9)  t  +  V*  {2  (f C  Ite)  t  cos  (n<t>)  +  2j(f<t>  l*9)  4>  sin  (n<{>) } 
~  u  n  n 


jt  =  (f ^  1^)  $  +  N  '  { 2 j  It^>)  t  sin  (n <t>)  +  2(f^  1^)  $  cos  (n<t>) . 
~  u  /  j  n  n 


where  ffc  and  f^1  are  row  vectors  of  the  f^(t)  and  fj(t),  respectively.  If  pf^(t) 

and  pfj(t)  are  the  triangle  functions  themselves,  rather  than  the  four-impulse  approx¬ 
imations  (A-31)  and  (A-32)  to  the  triangle  functions,  then 

oo 

I0i  *  +  ^2Ini  ^  COS  +  Jni  ^  Sin 


(A-55) 


l0i  *  + 


^  ^  {2jl^  t  sin  (n<j>)  +  21^  £  cos  (n4>) ) 


4.  NEAR  E-FIELD 

Near  g-fields  are  computed  in  a  manner  similar  to  that  described  by  Bevensee  [6]. 
s  t  h 

The  E-field,  E^,  radiated  by  the  n  mode  of  BOR  current  J  ,  along  a  circular  path 
concentric  with  the  BOR  axis  is  sought.  This  path  is  thickened  into  a  ribbon  of  sur¬ 
face  area  S  (Figure  A-8).  A  test  surface  current 
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-n 


~n 


(A-56) 


having  e  circumferential  variation,  resides  on  S  where  the  t  component,  J1" ,  lies 

in  the  plane  formed  by  the  BOR  generating  curve  and  BOR  axis  and  is  the  <(> -directed 

~n 

component.  By  reciprocity  [7],  the  t  and  <j>  components  of  E  (J  )  satisfy 

~n  -n 


ff  ES(J  )-J3ds=  //  E  (Ja)-J  ds  (A-57) 

J  J  '•n  '■n  ~n  J  J  ~ n  -n  -n 

S  s 


_  _ 

where  a  =  t  or  <J>  and  E  is  the  E-field  radiated  by  J  . 

~n  J  ~n 

With  given  by  the  nth  mode  terms  in  the  summation  (A-8) ,  the  right-hand  side 
of  (A-57)  becomes 
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(A-58) 
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where  I  and  I 
n  n 

t  h 

The  i  elements 


~a  t 

are  defined  following  (A-12)  and  Z 
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c  Zat  ,  «a<f> 
of  Z  and  Z  are 
n  n 


and  Za 
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are  matrix  row  vectors. 
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Figure  A-8.  Body  of  Revolution  with  Near-Field  Test  Ribbon 
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Figure  A-9.  Near-Field  Test  Ribbon  Triangle  Function  (Solid) 
and  Four-Impulse  Approximation  (Arrows) 
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where  a  =  t  or  4>  and  b  *  t  or  ♦  ,  and  the  vector  and  scalar  potential  functions,  A 
and  $,  are  given  by  (A-5)  and  (A-6 


Let  JC 

and  be  defined 

-n 

-n 

JZ  =  t 

£Ct)  e-jn* 

(A-60) 

j*  =  <»> 

'•n 

f(t) 

(A-61) 

where  pf(t)  is  the  "triangular"  function  shown  in  Figure  A-9.  The  right-hand  side  of 
(A-59)  then  is  within  a  minus  sign  of  (A-12)  if  m  and  n  in  (A-12)  are  each  replaced 
by  -n,  and  f^(t)  and  ff(t)  are  each  replaced  by  f(t).  In  this  comparison  of  (A-59) 
with  (A-12),  the  1 4>  and  <pt  terms  on  the  left-hand  side  of  (A-12)  are  reversed  from 
these  of  (A-15) ;  therefore,  the  four-impulse  approximations  to  pf(t)  and^-(pf(t)) 
are  given  by 
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where 
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and  comparison  with  (A-37)  through  (A-40)  and  (A-25)  through  (A-27)  yields 
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(A-63) 


(A- 64) 
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(A-66) 


In  (A-63)  to  (A-66) 
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=  j  d<j)'  c  kR  -  cos  (n<j> ' ) 


(A-67) 


Tl 

=  y  d<f>'  CRR  cos  <f>'  cos  (nd> ’ ) 


(A-68) 


*6  -  y  RR'  sin  <j>'  sin  (n«J>') 


(A-69) 


(P;L,  -  P  )2  +  (z^  -  z  )2  +  4p±p  sin2  (|-) 


(A-70) 


p '  =  p  +  4i  -  4 


i'  -  p  +  2i  -  2 


A-22 


v  is  the  angle  between  the  t  coordinate  of  the  test  surface,  S,  and  the  z  axis;  p 

-  -  q 

and  z  are  the  p,  z  coordinates  of  t  in  Figure  A-9. 

H  q 

Due  to  the  rotational  symmetry,  the  variation  of  J  (Equation  (A-9))  implie 

,  ind)  s  s 

the  same  e  variation  of  ?n(Jn)>  i*e.,  ^(J^)  may  be  written 


ES(J  )  =  tV  ejn*  +  ejn* 


~n  'n  n 


(A-71) 


where  E^  and  E^  are  <J>  independent.  If  Ea  is  approximately  uniform  over  the  width  of 
S,  the  left-hand  side  of  (A-57)  becomes  approximately 
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where  a  =  t  or  *.  Equations  (A-57),  (A-58) ,  and  (A-72)  yield 


(A-72) 


Ea  =  i  aat  ?c  +  za4>  p) 

n  4nA  n  n  n  n 


(A-73) 


The  t  and  <j>  components  of  the  E-field  E  radiated  by  the  BOR  surface  current 
are  expressed  as 


t-ES  -  E*  +  ^  (EC  ejn<,>  +  EC  e-Jn^) 

0  n  -n 


(A-74) 


♦  •  ES  =  E*  +  Vs  (E*  ejn*  +  E*  e“jn*) 
'  U  tm:  <d  n  -n 
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where  E^  and  are  given  by  (A-73).  The  positive  and  negative  "n"  modes  —  e.g., 

E^  and  EC  —  can  be  combined,  as  shown  below,  for  each  polarization,  0  and  $  ,  of 
incident  plane  wave  (A-l)  and  (A-2).  In  the  following  development,  a  second  super¬ 
script  on  E1",  E^,  IC,  or  1^  refers  to  this  incident  field  polarization, 
n’  n  n  n 

It  is  apparent  from  (A-63)  through  (A-69)  that 
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ZyT  = 

-n  n 


(A-75) 


From  (A-73),  (A-53),  and  (A-75),  therefore, 


Et9  -  Et6 
-n  n 


E*6  =  -E*6 


for  a  0  polarized  incident  field  and 


EC*  -  -EM 
-n  n 


E"  *  ET" 
-n  n 


(A-76) 


(A-77) 


for  a  <t>  polarized  incident  field.  Equations  (A-74),  (A-76),  and  (A-77)  yield 


Es  =  (e!;9  +  Ete  2  cos  n«>)t  +  T**  E*9  2j  sin  n<j>  $ 

0  1  ■*  n  £  n 


(A-78) 
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for  a  6  polarized  incident  field  and 

00  00 

Es  -  En^  2i  sin  n<t*  *  +  (Eq*  +  2  cos  n$)$ 

n=l  «-l 

for  a  <J>  polarized  incident  field. 
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5.  NEAR  H-FIELD 

Near  H-fields  (radiated  by  the  BOR  electric  surface  current)  are  computed  in  a 
manner  analogous  to  that  for  E-fields.  Again,  reciprocity  is  employed  in  the  deri¬ 
vation  and  the  test  current  (now  magnetic  rather  than  electric)  is  given  the  same 
"triangular"  t  variation  (A-62).  A  simpler  t  variation,  such  as  pulse,  could  have 
been  chosen  for  computing  the  H-fields  since  derivatives  of  the  test  current  then  do 
not  occur;  however,  the  choice  of  triangular  variation  maintains  a  consistency  through¬ 
out  the  overall  analysis.  The  H-field  derivation,  and  computer  programming,  then  can 
follow  closely  the  H-field  development  of  Mautz  and  Harrington  [1]. 

The  test  surface,  S,  of  Figure  A-8  now  contains  a  magnetic  surface  current,  M^, 
having  e  circumferential  variation.  By  reciprocity,  [7] 
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where  E  now  is  the  E-field  radiated  by  M  ,  J  is  the  nth  mode  surface  current  on  S, 
~n  ~  }  ~n'  ~n 

and  H  is  the  magnetic  field  radiated  by  J  .  From  [1  (Section  III)] 
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(r-p')  x  Mn(r')  ds 


(A-81) 


where  r  is  a  field  point  at  which  E^  is  evaluated  and  r'  is  a  source  point  (on  S)  at 

which  M  is  evaluated.  Primes  denote  source  coordinates  and  vectors.  Then, 

-n 

M  =  MC  t*  +  (A-82 ) 

~n  n  n 


and 


(A-83) 


(r-r’)  x  M  (r ' ) 

~  ~n  ~ 


(r-r’)  x  t’  M_ 


(t\<n  +  (r-r’)  x  tp'  hr  (t’,*') 


The  cross  products  on  the  right-hand  side  of  (A-83)  are  evaluated  by  expressing  all 
vectors  in  terms  of  unit  vectors  p,  <j>,  and  z  in  the  p,  <p,  and  z  directions, 
respectively. 


r  =  pp  +  zz  (A-84) 

r'  =  pp'  cos  ( <p'-<p )  +  4>p'  sin  (<£'-$)  +  zz'  (A-85) 

t'  =  p  sin  v'  cos  ($'-<!>)  +  <j>  sin  v'  sin  ( (J> *  — <f>)  +  z  cos  v'  (A-86) 

<p'  =  -p  sin  (<f> '  -  (fi)  +  <p  cos  ( <f> *  — d>)  (A-87) 

Equation  (A-85)  has  been  obtained  by  first  writing 

r'  =  p'p'  +  zz'  (A-88) 


and  then  using  Figure  A-6  to  express  p*  in  terms  of  p  and  <j>.  To  verify  (A-86),  use 
Figure  A-7  to  express  t'  in  terms  of  p'  and  z  and  then  use  Figure  A-6  to  express  p' 
in  terms  of  p  and  <j>. 

Substitution  of  (A-84)  through  (A-87)  into  (A-83)  yields 

(r-r')  x  M  (5')  =  (p  (-p'  cos  v'  +  (z'-z)  sin  v')  sin  (<j>'-4>) 

~  ~  ~n 

+  ip  (-(p-p'  cos  (<p'-<p))  cos  v'  -  (z'-z)  sin  v’  cos  (d> '  — d>) ) 

(A-89) 

+  zp  sin  v'  sin  ( <{> ' — 4>  )  }  M^(t',<t>')  +  (p(z'-z)  cos  ($'-<» 

+  ^  (  z'-z)  sin  (d> '  — <f>)  +  z(-p'  +  p  cos  (<J>'-<J>)) }  (t ' , d> ' ) 

n 

Let  v  be  the  angle  between  t  and  the  z  axis,  v  is  positive  if  t  points  away  from 
the  z  axis  and  negative  if  t  points  toward  the  z  axis.  Then 

t  «  p  sin  v  +  z  cos  v  (A-90) 


A-27 


and  the  t  and  <f>  components  of  (A-89)  become 


t*[(r-r')  x  M  (r * ) ]  =  — (p *  sin  v  cos  v'  -  p  sin  v'  cos  v 

~  ~n  - 

-  (z'-z)  sin  v  sin  v')  sin  M^(t',<t>')  -  {((p’-p)  cos  v  (A-91) 

-  (z'-z)  sin  v)  cos  (<J>'-<t>)  +  2p'  cos  v  sin2  (^  ^  ( t * , 4> 1 ) 

$*[(r-r')  x  Mn(r')]  =  {  ( (p  *  — p )  cos  v*  -  (z'-z)  sin  v')  cos  (4> '  — <>) 

(A-92) 

-  2p  cos  v'  sin2  (—  }  M^(t',4>’)  +  (z’-z)  sin  (<p’~<p)  M^(t',^') 

The  distance  |r-r'  |  appearing  in  (A-81)  is  the  square  root  of  the  sum  of  the 
squares  of  (z-z')  and  the  projection  of  (r-r')  in  the  xy  plane.  Hence, 


|r-r' !  =  y(z-z')2  +  p'2  +  p2  -2pp'  cos  (<|> '  — <f>) 

- - - -  (A-93) 

=  y( p-p' )2  +  (z-z' )2  +  4pp’  sin2  (^  2~~) 

An  integral  with  respect  to  <}> '  results  when  the  surface  integral  in  (A-81)  is 
iterated.  Because  this  integral  with  respect  to  d> '  is  an  integral  of  a  2ir  periodic 
function  of  4> '  over  the  period  2n,  <ji'  may  be  replaced  by  4> '  +  <j>  without  changing  the 
value  of  the  integral.  Substitution  of  (A-91)  through  (A-93)  into  (A-81)  yields 
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(A-95) 


R  =  y/ (p-p')2  +  (z-z')3  +  4pp'  sin^  Gy-) 


(A-96) 


All  the  iterated  integrals  in  (A-94)  converge  because  the  integrands  are  at  least  as 


well  behaved  as  R 


Mn(r')  =  l£(t',*')  or  l£(t'.*') 


I 


where 

MC  -  £  £(f)  e-Jn*' 
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M*  =  ♦  f(t')  e^"*' 
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and  pf  is  the  "triangle"  function  of  Figure  A-9 .  With  given  by  the  ntk 
in  the  summation  (A-8) ,  the  right-hand  side  of  (A-80)  becomes 
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-n  -n  ~n  n  n  n  n 


where  a  *  t  or  <p,  P  and  i’1*’  are  defined  following  (A-12),  and  Uat  and  Ua<* 
n  n  °  n  n 

row  vectors.  The  i1”^  elements  of  Uat  and  Ua^  are 

n  n 
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(A-97) 


mode  terms 
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are  matrix 
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where 
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(A-109) 


utj  =-(’-<■( (Pj ~ P±)  cos  -  (Zj-zi)  sin  vi>G2  +  kp^  cos  v^ 


UU  -  Jk(zj  -  2i)G3 


(A-110) 


v  is  the  angle  between  the  approximate  BOR  generating  curve  at  a°d  t^ie  7 

axis,  v  is  the  angle  between  the  t  coordinate  of  the  test  surface  and  the  z  axis, 
and  ,  z are  the  p,  z  coordinates  of  t^  in  Figure  A-9.  G^,  G 2>  and  G^  are  given 

by  (A-103) ,  (A-104),  and  (A-105)  in  which  G  is  given  by  (A-95)  with  R  of  (A-96) 
evaluated  at  (p,  z,  p' ,  z')  =  (p^  z.^,  p..  ,  z  ) . 

Due  to  the  rotational  symmetry,  the  ejn<^  variation  of  (equation  (A-9))  implies 

the  same  e-*n^  variation  of  H  (J  ):  i.e.,  H  (J  )  is  expressed  by 

'n  ~n  ~n  ~n 


H  (J  )  =  t  HC  e^11*  +  0  H*  ejnt| 


(A-lll) 


t  u) 

where  and  are  <p  independent. 

If  is  approximately  uniform  over  the  width  of  the  test  surface,  S,  and  is 
given  by  (A-97),  the  left-hand  side  of  (A-80)  becomes  approximately 


4A 


-  JJ Hn(Jn)'M®  ds  =  -H^  //  f(t)pd<)>dt  =  -H^  4itA 


0  0 


(A-112) 


where  a  =  t  or  <j>.  Equations  (A-80),  (A-98) ,  and  (A-112)  yield 


Ha  =  77  (U3t  +  U3'  ?*) 
n  4rrA  n  n  n  n 


(A-113) 


The  t  and  <p  components  of  the  H-field  H  radiated  by  the  BOR  surface  current  are 
expressed  as 
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(A-114) 


where  and  are  given  by  (A-113) .  As  with  ES  of  (A-74) ,  the  positive  and  nega¬ 
tive  "n"  modes  of  H  can  be  combined,  as  shown  below,  for  each  polarization,  0  and  <f>, 

t  d)  -+t  ~*<P 

of  incident  plane  wave  (A-l)  and  (A-2).  A  second  superscript  on  H^,  1^,  or  1^ 
in  the  following  development  refers  to  this  incident  field  polarization. 

It  is  apparent  from  (A-99)  through  (A-102)  that 

l3Ct  =  -15“ 

-n  n 


u ^  =  ut4> 

-n  n 


u**  - 

-n  n 


IT*  =  -lTH 

-n  n 


(A-115) 


From  (A-113),  (A-53) ,  and  (A-115),  therefore. 


HtS  =  -Hte 
-n  n 


(A-116) 


H*9  -  H*9 

-n  n 


for  a  0  polarized  incident  field  and 


Ht<fl  =  Ht(| 
-n  n 


(A-117) 


H4*  = 

-n  n 


for  a  $  incident  field.  Equations  (A-114),  (A-116),  and  (A-117)  yield 
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6.  APERTURE  COUPLING 


The  method  for  computing  the  near  scattered  fields  of  a  BOR  with  edges  (Sections 
1  through  5  of  this  appendix)  is,  in  theory,  directly  applicable  to  computing  the 
fields  penetrating  a  conducting  BOR  with  rotationally  symmetric  apertures;  e.g.. 
Figure  A-l.  The  internal  field  is  the  vector  addition  of  the  scattered  field  (Sec¬ 
tions  4  and  5)  and  the  impressed  field.  For  plane  wave  excitation,  the  impressed 
field  is  the  incident  field  E  given  by  either  (A-l),  (A-2),  or  a  linear  combination 
of  (A-l)  and  (A-2).  Each  BOR  model  component  of  E6,  at  a  field  point  of  interest, 
is  readily  computed  from  relations  similar  to  those  of  Section  2. 

As  often  occurs  with  well-shielded  bodies  (e.g.,  those  having  only  tiny  or  thin 
apertures),  the  scattered  and  impressed  fields  are  nearly  equal  and  opposite  at 
points  within  the  body.  The  vector  sum  of  these  fields  then  is  a  subtraction  of  al¬ 
most  equal  numbers  with  an  associated  loss  of  accuracy.  This  undesirable  computation 
is  avoided  by  invoking  an  equivalence  theorem,  due  to  Schelkunoff  [8],  whereby  the 
cavity  field  is  computed  as  a  radiation  problem  once  an  equivalent  aperture  excita¬ 
tion  is  determined.  This  theorem  states  that  the  original  problem,  containing  sources 
external  to  the  body,  can  be  replaced  with  an  equivalent  problem  having  only 
aperture  current  sources  for  excitation.  Figure  A- 10  illustrates  this  where  the 
original  problem  (Figure  10a)  is  equivalent  to  the  sum  of  the  two  problems  (Figures 
10b  and  c).  The  induced  surface  current  J  in  the  vicinity  of  the  aperture  is  found 
for  the  apertureless  (cover  aperture  with  conductor)  body  in  the  presence  of  the  ex¬ 
ternal  sources  (Figure  10b).  Then,  (-J  )  becomes  an  aperture  excitation  resulting 
in  the  external  fields  of  the  original  problem  (Figure  10c).  The  application  of  this 
theorem  to  the  moment  method  formulation  of  the  BOR  aperture  coupling  problem  has 
been  described  in  [9].  An  intermediary  aperture  tangential  E-field  computation  was 
employed  and  the  method  called  B0R3.  The  method  is  called  the  Equivalent  Aperture 
Excitation  Technique,  EAET,  in  this  report. 

A  modified  method  that  avoids  the  aperture  field  computation  has  been  developed 
[10];  however,  knowledge  of  the  aperture  tangential  E-field  when  the  aperture  is 
"empty"  is  useful  in  incorporating  the  Aperture  Loading  Technique,  ALT,  discussed  in 
the  main  body  of  this  report.  The  ALT  accounts  for  odd  shaped,  material  filled, 
thin  apertures  such  as  occurs  at  "sleeve"  type  joints  between  BOR  sections.  Thus, 
the  following  description  of  EAET  is  an  adaptation  of  [9]  to  the  development  in 
Sections  1  through  5  of  this  appendix. 
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The  apertures  are  removed  (covered  with  perfect  conductors).  The  induced  current 
on  this  "short-circuited"  BOR  is  expanded  as 


js-  y  (isc  jc  +  is*,  j*  ) 

nj  ~nj  nj  ~nj 


(A- 120) 


■*►  £  <+>  t  d) 

where  column  vectors  IS^  and  IS^,  composed  of  the  coefficients  IS^  and  IS^  ,  respec¬ 
tively,  are  the  solutions  1^  and  I*  of  (A-20). 

The  apertures  are  replaced  (shorting  conductors  removed)  and  are  filled  with  the 
equivalent  aperture  current  source,  JA,  expanded  as 


JA  -  V*  (IAC  JC,  +  IA*  J*.  ) 
L-!  np  ~njn  np  ~ni 


(A-121) 


where 


IA*  -  -  IS* 
np  nj 


IAy  -  -  IS* 
np  nj 


(A-122) 


and  the  p  t-directed  (or  (^-directed) aperture  expansion  function  is  the  j  corre- 
sponding  generating  curve  expansion  function.  The  external  incident  field,  Ee,  is 
removed  and  this  new  problem  analyzed  as  follows.  First,  the  aperture  tangential 
component  of  E-field  is  expanded  as 


EA 

-  tan 


=  V'  (EA11  t  +  EA*  $) 
/  ,  np  np 


(A-123) 


where 


t  -  ejn* 

VAC  <$(t-t„ .  .  )— 5 — 

np  2jp+l  2it 


VAnp  a<«ij  +1>V 

P 


(A-124) 


6(t)  (Section  1  of  this  appendix)  denotes  the  unit  impulse  function,  and  t~  +1 

.  Z3P 

(Figure  A-3)  is  the  BOR  generating  curve  t  coordinate  of  the  triangle  function 
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pfj  (t)  (or  pf?  (t))peak.  The  generating  curve  triangle  function  is  the  p*"*1 
Jp  JP  P 

aperture  located  triangle  function. 

Each  VA1"  and  VA^  can  be  thought  of  as  port  voltage  sources  through  which  a  sur- 
face  current  flows.  One  seeks  values  of  VA^  an':’  VA^  for  which  J,A,  given  by  (A-121) 
and  (A- 122),  flows  through  the  aperture  ports.  These  values  are  obtained  by  relating 
the  IA1"  and  IA^  of  (A-121)  to  the  VAC  and  VA^1  by  the  matrix  equation 


np 


np 


np 


np 


a*  1 

n 

II 

U 

< 

1 _ 

[YA^]' 

VAC 

n 

-*■  p 
IAV 
n 

i 

rt 

(YA^] 

_ 1 

(A-125) 


and  solving  (A-125).  In  (A-125),  the  p  elements  of  IA^,  IA^,  VA^,  and  VA^  are 

IA*"  ,  IA^  ,  VAC  ,  and  VA^  ,  respectively.  The  matrix  of  [YA3^],  ab  =  tt,  <f>t,  1 4> , 
np ’  np ’  np  *  np  r  3  n 

and  (f>4>  submatrices  is  called  the  "aperture  admittance  matrix." 

f"  Tn  jjk 

The  pq  elements  of  [YA3  ] ,  ab  =  tt,  $t,  tip,  and  pp ,  are  found  by  "exciting" 

the  BOR  with  the  tangential  aperture  field,  EA  ,  given  by  (A-123).  This  requires 
satisying  the  boundary  condition 


EA 

~  tan 


on  S 


(A-126) 


g 

where  (Section  1  of  this  appendix)  E  is  the  field  radiated  by  the  BOR  surface  cur¬ 
rent,  J.  A  development  similar  to  that  following  (A-3)  applied  to  (A-126)  yields 
(A-20)  with  the  exception  that  the  jt  element  of  and  are,  respectively. 


VA 


np 


VA* 


np 


j  =  jP 
otherwise 

otherwise 


(A-127) 
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The  elements  of  1^  and  in  (A-20)  corresponding  to  aperture  expansion  functions  are 

t  rf) 

the  elements  of  IA  and  IA  ;  l.e.* 
n  n 


I  -  IA 
nj  np 


if  j  =  j. 


(A-128) 


r.  -  ia^ 

nj  np 


if  j  -  j. 


In  light  of  (A- 125),  (A-127),  and  (A-128)  and  defining  the  inverse  of  the  impedance 
matrix  in  (A-20)  as 


[y“] 

[Y**l 

•u 

N) 

1 _ 

[z'*l  ' 

n 

n 

_ 

n 

n 

w  _ 

[Z*c] 

n 

[ZH] 

n 

the  pqtl*  elements  of  [YA^*"],  [YA^C],  [YA^],  and  [YA^]  become 


itt:  = 

YCt  . 

npq 

nJPJq 

= 

Y*J  . 

npq 

njpJq 

» 

Yfc*  , 

npq 

nj  j 

JpJq 

npq 

(A-130) 


The  solution  of  (A-125),with  IA^  and  IA^  given  by  (A-122) ,  provides  the  aperture 
port  voltage  sources  that  correspond  to  the  aperture  tangential  E-field  within  the 
approximation  (A-124).  Excitation  of  the  BOR  with  these  voltage  sources  yields  the 
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